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Surface  and  edge  energy  of  electron  gas  in  nanocrystals. 

M.  V Entin  and  M.  M.  Mahmoodian 

Institute  of  Semiconductor  Physics,  Siberian  Branch  of  the  RAS, 
630090  Novosibirsk,  Russia 


Abstract.  The  energy  of  a many-electrons  nanoparticle  depends  on  its  form  via  the  surface  and 
edge  contributions  to  the  thermodynamic  functions.  The  surface  energy  of  free  electron  gas  is 
found  for  a semiconductor  with  the  Dirac  spectrum.  The  edge  energy  of  electron  gas  is  calculated 
for  quadratic  energy  spectrum  of  electrons  without  and  with  anisotropy.  We  find  also  spin-orbit 
corrections  to  the  surface  energy,  which  behave  like  the  inverse  square  of  crystal  size. 


Quantum  dots  with  the  large  number  of  electrons  are  intermediate  objects  between  the 
quantum  and  the  classical  limits.  On  the  one  hand,  their  sizes  are  large  compared  with  the 
electron  wavelength.  On  the  other  hand,  the  transport  and  thermodynamic  properties  of 
the  system  are  affected  by  the  quantization  of  the  electron  states. 

In  the  degenerated  electron  system  the  surface  contribution  to  the  thermodynamics  of 
small  particles  is  determined  by  the  small  ratio  of  the  Fermi  wavelength  of  the  electron 
1 /kf  to  the  particle  size  L.  The  presence  of  the  surface  leads  to  the  regular  corrections  in 
powers  of  this  parameter  to  the  chemical  potential  of  electrons  [1]. 

The  surface  contributions  to  the  energy  of  the  electron  gas  lead  to  the  different  physical 
effects.  In  particular,  they  affect  the  surface  tension  in  small  particles  and,  consequently, 
their  equilibrium  form.  Establishment  of  the  equilibrium  between  the  electron  gas  in 
microparticles  with  different  size  or  form  is  accompanied  by  their  spontaneous  charging  [ 1 ]. 
As  was  shown  in  [2],  owing  to  charge  discreteness,  the  chemical  potentials  in  microparticles 
are  not  completely  equalized.  As  a result,  the  system  of  metal  granules  becomes  a gapless 
insulator  (the  gapless  Hubbard  insulator). 

Aside  from  the  surface  contribution  to  the  energy,  in  facetted  nanocrystals  there  are 
contributions  caused  by  their  edges  and  vertexes  [3].  The  goal  of  the  present  report  is  the 
study  of  surface  and  edge  energies  of  three-  and  two-dimensional  facetted  nanocrystals. 

Edge  contribution  to  the  energy  in  the  free  electron  model 

Here  we  consider  two  or  three-dimensional  gas  of  free  electrons,  bounded  by  an  angle 
0 < ip  < 0,  r — Jx2  + y2  < R (the  dihedral  angle  in  3D  case).  The  electron  wave 
function  \j/  satisfies  the  condition  \[r  — 0 on  the  angle  borders.  We  search  the  corrections 
to  the  Q potential  at  zero  temperature,  caused  by  presence  of  the  edge  and  the  boundary. 
The  contribution  to  the  potential  from  the  vicinity  r < ro  of  the  edge  r — 0 is 

Lz  r0  <t> 

& = - j dz  j rdr  J dcp^j^-  - k^/2me)\jrnym(z,  r , <p)  |2  (1) 

0 0 0 

Here  sn_m  is  the  energy  level  of  the  electron  with  the  main  n and  the  magnetic  m quant  um 
numbers,  kz  is  the  momentum  along  the  edge,  me  is  the  electron  effective  mass,  fi  — 
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k2F/2me  is  the  chemical  potential.  The  axis  z is  directed  along  the  edge  and  Lz  is  the  edge 
length.  In  2D  case  this  variable  and  corresponding  integration  should  be  excluded. 

When  the  radius  R of  the  sector  goes  to  infinity,  £2  ceases  to  depend  on  the  sector  size. 
In  the  limit  kpro  1 £2  yields  an  asymptotic  expansion  (3D  case): 


•y 

£2  = cojLz4>r0/2  + a>2Lzro  + w\Lz  + o(l). 
The  wave  function,  satisfying  boundary  conditions  is 


0 n,m(z,r , (p) 


exp (ikzz)Jv(kr)  sin(ixp). 


Tcm 


0 


(2) 


(3) 


Using  the  wave  function  (3)  we  find 

kF  l oo 

£2  = - — [ dk-rU 4 f xdx[  lnx  H — (1  — v2)]  jlm(ux).  (4) 

merk  JT  J J 2 , * 

0 0 0 m=1 

Here  u = (kj  - k2)'/2r(). 

The  expansion  of  (4)  in  the  orders  of  ro  at  ro  — >•  oo  gives  o>3,  0)2  and  co\.  The  term 
0)3  = —k5f/\5TC2me  is  the  usual  £2  potential  of  3D  gas.  The  term  a>2  = kF/327Tine 
represents  the  surface  contribution  [1,2]. 

The  edge  contribution  a)\  yields 


CO]  (4>)  = 


3 6 JT  /?7 1,  v 0 


(5) 


Physically  the  dependence  of  co\  (0)  is  explained  by  the  same  factors  as  the  positive  sign 
of  the  surface  contribution  to  the  energy.  Zero  boundary  conditions  on  the  crystal  surface 
deplete  the  electron  gas  density  in  the  boundary  region.  As  a result,  the  size  of  occupied 
region  decreases  that  raises  the  energy  of  electron  gas  at  the  unchangeable  number  of 
electrons.  The  volume  of  the  depletion  region  grows  with  the  angle  that  leads  to  the  same 
change  of  the  energy  of  electron  gas  if  both  crystal  volume  and  surface  area  don’t  vary. 
The  change  of  co\  (0)  sign  at  0 — tt  is  caused  by  extracting  of  the  surface  part  of  energy. 

Two  cases  are  of  special  interest.  They  are  the  right  angle,  typical  for  the  facetting  of  a 
cubic  crystal  and  the  angle  cp  = 2n,  corresponding  to  the  rift  in  the  crystal.  In  these  cases 
o>i  (tt/2)  = —kp/24nme  and  co\  (2tt)  = k2F/24june. 

By  means  of  (5  ) the  number  of  electrons  in  nanocrystal  of  volume  V . the  surface  area 
S and  edges  with  angles  0„  and  lengths  Ln  is  expressed  through  the  chemical  potential: 


3 i 

3£2  (2mefJ.)2V  me[xS  ( 2me[i )2  ^ n 0n  \ 

dii  3 tt2  4tt  \2n  “ " v0„  Tt  / 


(6) 


The  analogous  formula  is  fo  und  for  the  edge  energy  of  a 2D  quant  um  dot  with  a polyg- 
onal border.  The  angle  energy  of  2D  gas  is 


ft)l  (0)  - - 


kF  (n_  _ 0\ 
A d>  7 r / 


24m e V 0 jt 


(7) 
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Anisotropic  electron  energy  spectrum 

The  expression  for  the  edge  contribution  to  the  energy  allows  generalization  to  the  case  of  the 
anisotropic  quadratic-law  energy  spectrum  s(k)  — J2  The  affine  transformation 

x-  — Xi(mj/me)1/2,  where  me  — converts  the  Schrodinger  equation  to 

isotropic  one,  for  which  the  previous  consideration  is  valid. 

If  we  establish  the  dihedral  angle  by  means  of  the  normal  vectors  ni  and  n2  to  the  forming 
planes,  the  transformed  angle  between  these  planes,  entering  into  previous  formulas,  is 
defined  by  the  relation: 


3 

cos  4>  = y 

i= 1 


«W»2i 

nij 


i 

7 


and  the  edge  length  Ln 


r 3 -.1/2 

2_,  Llni  is  changed  on 

L ;=1  -1 


nij  jr  2 
me  ^ ni 


(8) 


Dirac  model 

The  Dirac  model  is  the  simplest  Hamiltonian,  describing  the  energy  spectrum  of  a semicon- 
ductor. The  most  adequately  it  corresponds  to  a cubic  semiconductor  with  a non-degenerate 
central  main  minimum  of  the  conductivity  band  and  maximum  of  the  valence  band.  The 
standard  Lattinger  model  with  degenerate  maximum  of  the  valence  band  turns  into  the 
Dirac  model,  if  an  anisotropy  of  the  band,  heavy  holes  and  spin-orbitally  splitted  band  are 
overlooked.  In  comparison  with  the  quadratic-law  spectrum  the  Dirac  model  allows  to  take 
into  account  non-parabolicity  of  the  electron  energy  spectrum. 

The  Dirac  equation  for  the  spinor  component  u and  v has  the  form 

(s  — mec2)u  — c(ka)v  = 0,  (s  + tnec2)v  — c(kcr)u  — 0.  (9) 

Here  s is  the  energy  of  the  electron,  counted  out  from  the  center  of  the  forbidden  band  with 
the  width  Eg  — 2 mec2,  a is  the  Pauli  spin  matrix. 

Each  of  the  spinor  component  satisfy  the  Klein-Gordon  equation 

(s2  — ( mec2)2)u  — c2k2u  — 0.  (10) 

The  equation  (10)  transforms  into  the  Schrodinger  equation  after  substitution  (e2  — 
(mec2)2)  —*■  2 mec2E.  If  we  suppose  zero  boundary  conditions  for  the  large  components 
of  the  wave  function  u — 0,  the  agreement  with  the  Schrodinger  equation  becomes  full. 
Hence,  the  previous  formulae  prove  correct  if  we  change  /x  = k2F/2me  —*■  y+EEs  referring 
ix  to  the  conduction  band  bottom. 

Spin-orbit  induced  surface  corrections 

The  spin-orbit  interaction  is  relativistic  and  usually  small  effect.  Nevertheless,  spin-orbit 
interaction  is  known  [4-5]  to  be  amplified  due  to  confinement  by  quantum  wells.  If  the 
width  of  a well  tends  to  infinity,  while  the  Fermi  energy  keeps  constant,  the  corrections  to 
the  electron  spectrum  go  down  inversely  as  the  size.  This  give  rise  a new  mechanism  of 
the  surface  corrections  to  the  thermodynamic  functions. 
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Let  us  consider  a planar  sample  with  the  width  L and  the  normal  z.  The  spin-orbital 
correction  to  the  Hamiltonian  yields 


k 1 

Hso  = — — 
meL 


(cxxky  - cfykx)  X = x 


4 m2c2 


(ID 


where  x is  some  constant  of  order  of  unit,  depending  on  the  specific  surfaces.  If  the  surfaces 
are  equivalent  y vanishes,  but  in  the  general  case  x is  finite. 

The  Hamiltonian  (11)  has  the  energy  spectrum  e(k)  = k2 /2me  ± j^Jk 2 +k2. 

Substituting  this  spectrum  to  Q potential  we  find  the  correction  to  the  number  of  electrons 
in  the  crystal  AN$o  = V'/2L  2k5l./5jt2.  The  found  correction  has,  evidently,  surface 
origin.  So  the  order  of  such  correction  will  be  the  same  in  a nanocrystal  of  any  shape  with 
the  characteristic  size  L.  This  correction  turns  out  to  be  weaker  than  the  surface  correction, 
proportional  to  L-1,  but  may  be  comparable  with  the  edge  corrections. 


The  numerical  calculations  of  the  energy  states  in  many-electron  quantum  dots  are 
complicated  enough  problem  because  of  the  exponential  increasing  ofthe  number  of  energy 
states,  included  in  the  Hamiltonian  matrix,  with  the  growing  number  of  electrons.  At  the 
same  time  the  considered  approach,  based  on  the  expanding  of  thermodynamic  quantities 
in  powers  of  size,  allows  to  get  simple  estimations  without  complex  calculations. 
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